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Abstract 

Within the context of finite deformation elasticity theory the prob- 
lem of deforming an open sector of a thick-walled circular cylindrical 
tube into a complete circular cylindrical tube is analyzed. The analysis 
provides a means of estimating the radial and circumferential residual 
stress present in an intact tube, which is a problem of particular con- 
cern in dealing with the mechanical response of arteries. The initial 
sector is assumed to be unstressed and the stress distribution result- 
ing from the closure of the sector is then calculated in the absence of 
loads on the cylindrical surfaces. Conditions on the form of the elas- 
tic strain-energy function required for existence and uniqueness of the 
deformed configuration are then examined. Finally, stability of the re- 
sulting finite deformation is analyzed using the theory of incremental 
deformations superimposed on the finite deformation, implemented in 
terms of the Stroh formulation. The main results are that convex- 
ity of the strain energy as a function of a certain deformation vari- 
able ensures existence and uniqueness of the residually-stressed intact 
tube, and that bifurcation can occur in the closing of thick, widely 
opened sectors, depending on the values of geometrical and physical 
parameters. The results are illustrated for particular choices of these 
parameters, based on data available in the biomechanics literature. 

Keywords: nonlinear elasticity; finite elasticity; residual stress; elastic tube; 
existence and uniqueness; stability 



1 Introduction 



Residual stresses have a very important role in the mechanical functioning of 
arteries and the existence of residual stresses is well documented, as in, for 



example, Chuong and Fung (1986), Vaishnav and Vossoughi (1983) and Fung 



(1993); see also the review by Humphrey (2002). They are demonstrated by 



the simple device of cutting radially a short ring of an artery, the result of the 
cut being the springing open of the ring into a sector, thereby releasing some 
residual stresses; in general, some residual stress remains, however, as shown 
by Vossoughi et al. (1993) and Greenwald et al. (1997); see also the review 



by Rachev and Greenwald (2003). A crude estimate of the circumferential 
residual stresses is obtained by measuring the resulting angle of the sector 
into which the ring deforms, although it has to be said that reliable quan- 
titative data remain elusive. This so-called opening angle method has been 
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analyzed in some detail by several authors, including, for example, 


Delfino 


et al 


• ( 


L997) 


Zidi et al. (1999), Rachev and Hayashi 


'1999) 


Holzapfel et 


al. (; 


>00 


0), 


gden and Schulze-Bauer (2000), Matsumoto and Sato 


(2002). 


Ogden 


2003 


), Raghavan et al. (2004) and Olsson et al. 


(2006 


. The typical 



approach is to assume that the (unloaded) opened sector is free of stress 
and circular cylindrical and then to construct a deformation that brings the 
sector into an intact tube and to calculate the resulting radial and circum- 
ferential (residual) stress distributions. This has been done for a single layer 
and for two-layered tubes. In general this method does not account for any 
stress or deformation in the axial direction. However, in reality, not only 
does the ring open into a sector, but the length of the arterial segment tends 
to change and there are also residual stresses in the axial direction, although 
quantitative information about axial residual stresses is distinctly lacking. A 



recent 3D analysis by Holzapfel and Ogden (2010) is the first attempt to cal- 



culate the combination of radial, circumferential and axial residual stresses. 
The purpose of the present paper is to provide an analysis of the opening 
angle method, with particular reference to the questions of existence and 
uniqueness of the residually-stressed configuration and its stability. 

In Section [2] we review the description of the geometry of the deformation 
which first takes the sector into an intact cylinder, allowing for a possible 
length change, and then allows inflation under internal pressure combined 
with axial load while maintaining circular cylindrical symmetry. The corre- 
sponding stress components and equilibrium equations are then summarized 
in respect of an isotropic elastic material whose properties are described in 
terms of a strain-energy function. The existence and uniqueness of the re- 
sulting configuration is then examined in Section [3j under the assumption 
that the strain energy is convex as a function of a suitably chosen deforma- 
tion variable. The Mooney-Rivlin and the Fung strain energies obey this 
assumption, and we use them, as well as some experimental data on arteries, 
to illustrate the analysis. Section [4] is devoted to an analysis of the stability 
of the residually stressed tube (in the absence of lateral loads) on the basis of 
the theory of small incremental deformations superimposed on a finite defor- 
mation, and for this purpose an appropriate version of the Stroh formulation 
is adopted. For simplicity of illustration, attention is restricted to prismatic 
incremental deformations. We then treat numerically the case of tubes made 
of solids with the Mooney-Rivlin strain energy, and compare the predictions 
to some simple experiments we performed with silicone rubber. 



3 



2 Problem formulation 



2.1 Geometry of the deformation 

We consider an annular sector of a circular cylindrical tube with geometry 
defined by 

A < R < B, -(2tt - a)/2 < < +(2vr - a)/2, < Z < L, (1) 

where A and B are the radii of the inner and outer faces, respectively, of the 
sector, L is its length, and (R, 9, Z) denote the cylindrical coordinates of 
a point in the material, with corresponding orthonormal basis vectors (En, 
E®, E z ). Here a G (0,27r) is the so-called opening angle. The sector is 
assumed to be unstressed in this reference configuration, which we denote by 

The sector is then deformed into an intact (circular cylindrical) tube so 
that the plane faces originally at = ±(7r — a/2) are joined perfectly to- 
gether and there is an accompanying uniform axial stretch A 2 = l/L > 1, 
where I is the deformed length of the tube. We refer to this new configura- 
tion as the residually-stressed configuration, which we denote by B r . In this 
configuration there is no traction on the curved surfaces of the tube, but in 
general axial loads are required to maintain the deformation since it turns 
out that the axial stress depends on the radius r and cannot therefore be 
set to zero pointwise on the ends of the tube. Finally, the tube is subjected 
to a uniform internal pressure of magnitude P per unit deformed area and 
is maintained in a new circular cylindrical configuration, which we refer to 
as the loaded configuration, denoted B\, the three separate configurations are 
depicted in Figure [1} 

Let (r, 9, z) denote cylindrical coordinates in either the residually-stressed 
configuration B r or the pressurized configuration Bi, with corresponding or- 
thonormal basis vectors (e r , eg, e z ). The deformation may be described by 
the equations 

r = r{R), 9 = kB, z = \ z Z, (2) 

from which we may calculate the deformation gradient tensor F = Grad;c 
as 

F = r'e r ®E R + (kr/R)e e ®E e + \ z e z ® E z , (3) 
where the prime denotes differentiation and 

2tt 



Figure 1: Depiction of the three configurations: (a) unstressed open sector 
(configuration B ), with the opening angle a shown as n/2; (b) residually 
stressed and stretched tube (configuration B r ); (c) tube in (b) subject to 
internal pressure P (configuration Bj). 

is a measure of the opening angle. The left Cauchy-Green tensor is then 
calculated as 

B = FF T = r' 2 e r <g> e r + {kr/R) 2 e e ® e e + X 2 z e z ® e z , (5) 

and it follows immediately that the principal axes of B are e r ,e e , and e z . 
The corresponding principal stretches are 

Ai = r', A 2 = %■, A 3 = \ z . (6) 
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Arterial wall tissue is generally considered to be essentially incompressible, 
so that the constraint det F = 1 is enforced. For the present deformation 
this constraint yields 

where a = r(A) is the radius of the inner surface of the tube in configuration 
B\. Similarly, we define b = r(B) for the outer surface. It then follows from 
from ^ that the principal stretches become 

Notice from ([7]) that if a 2 k\ z = A 2 then the principal stretches are con- 
stants. In other words, the deformation is homogeneous, of the form 

r = R/y/kK, 9 = kQ, z = \ Z Z, (9) 

with principal stretches A x = 1/ \/k\ z , A 2 = \/k/X z , A 3 = A z . 

2.2 Stress components and equilibrium 

Henceforth we restrict the analysis to homogeneous, incompressible, isotropic 
hyperelastic solids. The strain-energy function is then a symmetric function 
of the principal stretches, which we write as W(\i, A2, A3), subject to the 
constraint A1A2A3 = 1. For such materials the Cauchy stress tensor cr is 
coaxial with B and for the considered deformation the principal Cauchy 
stress components are 

BW BW , dW , , 

a rr = -q + Ai— -, a ee = -q + X 2 ^—, a zz = -q + A 3 — -, (10) 
(Mi aA 2 aA 3 

where q is a Lagrange multiplier introduced by the constraint of incompress- 
ibility, and Ai, A 2 , A 3 are given by rt8J). 

Because the stretches are independent of 9, and z the 9 and z components 
of the equilibrium equation ensure that q is also independent of 9 and z. The 
only non-trivial equation remaining is the radial equation 

— a rr + -(cr rr — oee) = 0. (11) 
dr r 



6 



With this equation we associate the boundary conditions 



a rr (a) = -P, a„{b) = 0, (12) 

where P is the internal pressure in configuration Bi and the outer boundary 
is free of traction. 

At this point it is convenient to introduce the notations 

a 2 r 2 b 2 B 2 

x a = k\ z — , x = k\ z — , x b = k\ z — , e = — - 1 > 0, (13) 

e being a measure of the thickness of the initial annular sector. Then it 
follows easily that 

" = ^T- < 14 > 

For any fixed value of the axial stretch X z we may now regard W as a 
function of x through the stretches, which are given in terms of x by 

A -vfc' A -/f' (15) 

together with A 3 = X z , and we write W(x) = W(l/ y/kX z x, ^kx/X z , X z ). 
From ( 10 ) it then follows that 

dW dW - 
a ee - a rr = X 2 — - Ai— = 2xW (x) . (16) 
oX 2 oXi 

Noting that 

=2x(l- x), (17) 

dr 



we may integrate equation (11) and apply the boundary conditions Q12j) , and 
thus conclude that ^ 

fx b \y' (x) 

P= / ^-^dx. (18) 
.L I- x 



This equation, together with (14), enables the location of the deformed inner 
radius to be determined for given values of the pressure P, the thickness (as 
measured by e), the opening angle (as measured via k), and the axial stretch 
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X z . For completeness we compute the corresponding stress distribution, given 
by 



W (s) dW dW 

ds, a ee = a„ + 2xW (x) , a zz = a rr + A 3 — Ai— — . 

1 — s 0X3 0X1 

(19) 

In the special case of the homogeneous deformation (p)]), we have x = 1 
and we find from (16) that 



a rr = 2W'(1) = constant, 



so that we can integrate the equation of equilibrium (11) to give 

2W'(l)\n (r/b), 



(20) 



(21) 



where the boundary condition (12 (2 has been used. The boundary condition 
( 12 )i then yields 



P = W\l) In (b/a) = W'{1) In (B/A) . 



(22) 



This is the pressure that needs to be applied inside the tube in order to 
produce a homogeneous deformation through the tube wall. We note that in 
this case, both a rr , given by (21), and aee, given by 



P[l + ln(r/6)]/ln(6/a) 



(23) 



have logarithmic variations. Hence, in the case of homogeneous deformation, 
the stress is a slowly varying function of the radial coordinate. 

It is clear that in general a zz depends on r. In the case of homogeneous 
deformations this follows immediately from (21) and (19)3. Thus, we em- 
phasize that the pointwise end condition a zz = cannot be adopted and the 
residually stressed configuration B r is not load free. In some situations it may 
be possible to enforce the condition that the resultant axial load vanishes, 
but here we merely consider that the axial stretch is fixed, so that the axial 
load has to be adjusted accordingly to accommodate this. 



2.3 Convexity of the strain-energy function 

Because we are considering an incompressible isotropic elastic material, the 
strain-energy function can also be written as W = W (Ji, I2), where I±, I2 are 
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the first and second principal invariants of the Cauchy-Green deformation 
tensor B. Expressed in terms of the principal stretches these are 



h — Aj + A 2 + A3, 1% — \\\\ + X\\\ + A2A3 
From (J8I) we then have explicitly 
R 



(24) 



h 
h 



k\ 7 r 



2 / ; \ 2 
/ kr 



+ 



1 (R 

A| + Vkr 



\R 



+ 



+ A* 



1 kcc 2 

T\ 1" "a — I" A~, 

/cA,x A, 



fcA 2 r 



^ + k\ z x + A 2 



-2 



(25) 



The function W (x), occurring in the integrand of (18), can therefore be 
written as 



W (x) 



k 2 x 2 -1 dW ^ 2 dW 
+ A' 



k\?x 2 



dh 



dh 



(26) 



Assuming that the standard empirical inequalities dW/dli > 0, dW/dI 2 > 
hold, an immediate conclusion is that independently of the form of W (Ji, I2) 



W' (x) — a; = x r 



1/k = (2tt - a)/(27r). 



(27) 



We assume here and henceforth that W (x) is a strictly convex function, 
i.e. W" (x) > 0. It follows that for any given value of X z , x = x m = 1/k 
yields the unique minimum of the strain energy function. Further, by (24) 
we have 

h(x m ) = 2A; 1 + Xl I 2 {x m ) = 2\ z + A; 2 , (28) 

and the minimum value of the strain energy is therefore independent of the 
measure of the opening angle k = 2ii/ (2tt — a). 
From equations (11), (16), and (17) we have 

W'(x) 



do~ rr 

dx 



1 — x 



(29) 



and consequently, the radial stress a rr also has a unique extremum, a min- 
imum, at x = x m . The corresponding relation for the hoop stress aee is 
derivable from ( 19 ) and has the form 

dam W'(x 



dx 



+ 2W'(x) + 2xW"[ 



x 



X 



(30) 



9 



In contrast to the radial stress, the determination of its singular point (s) 
depends on the form of the strain energy function. Clearly, the hoop stress is 
increasing at x — x m , the radius where the radial stress is at a minimum. At 



that point, according to (19) 2 and (27), o~oq = a rr and, by l\8ti, Ai = A 2 ; thus 



the in-plane principal stresses and the stretches are equal there. We refer to 
the surface defined by the neutral surface. The associated radius 

(r or R) depends, in general, on both A^ and k. 



Another consequence of (26) and convexity is that 



W (x) < for < x < x m = l/k, 
W (x) > for x > x m = l/k. 
Application of these conditions to the tube, however, depends on whether 



x a > Xb or Xb > x a . It is easy to show from (14) that either (i) x a > x& > 1 
or (ii) 1 > Xb > x a . In case (i) we have x m = l/k < 1 so x m does not lie 
within the required range and there is no neutral surface in this case. In case 



(ii) we require Xb> x m > x a . Thus, (31) can be made more explicit 



W (x) < for x a < x < x m = l/k, ^ 
W (x) > for 1 > x b > x > x m ee l/k. 

This will be used in the next section to determine existence and uniqueness 
of the unloaded, residually-stressed configuration. 

In closing this section, we note that convexity of W(x) is a common fea- 
ture of many standard strain-energy functions, including the Mooney-Rivlin 
model, 

W M R= ( ^(h-S) + I ^(I 2 -3), C>0, D>0, (33) 

where the material constants C and D have dimensions of stress, the Gent 
model 

fJ-Jm , /, h — 3 



W G = -^\n[l-^—), M >0, J m >0, (34) 

*J m. 



and the Fung model 



W F = ^[e c ^-l], M >0, c>0, (35) 

where /i is the shear modulus in the unstressed configuration and J m and c 
are dimensionless constants. Convexity may easily be checked. For instance, 
Wj[J R (i) = (C + DX 2 z )/(kX z x 3 ) which is clearly positive. Similar, but longer, 
expressions can be found for Wq and Wp, which are also strictly positive. 
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3 Existence and uniqueness of the residually- 
stressed configuration 

Here and henceforth, we study the residually-stressed tube in configuration 
B r . Setting P = in pi yields 







j' b 



W'(x) 



dx, 



x 



£ + Xg 

e + 1 ' 



(36) 



This equation determines an inner radius of the unloaded residually-stressed 
tube, measured by x a , for a given measure k of the opening angle and given 
axial stretch X z . 

Two questions need to be answered. 



i) Does equation (36) have a positive solution x a ? In other words, is it 



always possible to close a cylindrical sector in the way indicated when 
no load is applied to its curved faces? 

(ii) If a solution exists, is it unique? 

First consider the possible existence of a homogeneous residually-stressed 
configuration, when the deformation is in the form ([9]), which is equivalent to 
assuming that x a = 1. Then, by (22), vanishing of P is equivalent to W(l) = 
0. According to (27) this leads to k = 1, which is not admissible since k > 1. 



Thus, there cannot be a homogeneous residually-stressed state of deformation, 
i.e. a residually-stressed configuration is necessarily inhomogeneous. This, 



of course, is well known in general, as shown by, for example, Hoger (1985); 



see also the discussion in Ogden (2003) 



Assume now that x a > 1. Then, as indicated following (31), we must 
have x a > Xb > 1 > x m , and hence, according to (31), the integrand in (36) 



has the same sign over the entire range of integration, and there is therefore 



no solution to (36) for x a > 1. 



For a solution x a of (36) to exist, it must lie in the range < x a < 1. If 
either x a < Xb < x m or x m < x a < Xb then there is again no solution to (36), 
because the integrand is one-signed over the range of integration. Thus, we 
must have the ordering x a < x m < Xf> and (32) is applicable. Expressing x& 



in terms of x a by means of (36)2 in this latter inequality yields the following 
range of values for x a : 



[1 - (k - l)e]x m < x a < x r , 



(37) 
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Then, to ensure that the solution x a to (36) is positive, we enforce the in- 
equality e < l/(k — 1), i.e. 



B 2 



< 



k 



A 2 k 



2vr 



(38) 



which prescribes a maximum ratio of outer to inner radius for the original 
annular sector for a given opening angle, or equivalently a maximum value 
of the opening angle for a given radius ratio. 

Now, addressing question (i) above is equivalent to determining whether 
the function / defined over the range (37) by 

«(e-n/)/(e+i) w'(x 



f(y) 



1 — X 



dx, [1- (k- l)e)x m < y < x T 



(39) 



has a zero. First we examine the value of f(y) at the upper bound of the 
range, which is 

"(e+x m )/(e+l) W'( x \ 



f(Xm) 

and this is positive because 



-dx, 



x 



e + 1 



(40) 



(41) 



so that the integrand in (40) is always positive over the range of integration. 
Next, we examine f(y) at the lower bound of the range (37); after some 
simplification, we find that 



f([l-(k-l)s]x r , 



[l-(fe-l)e]a:„ 



W'{x) 

1 — X 



dx, 



(42) 



which is strictly negative because 

[1 - (k - l)e}x m <x m <l, 



(43) 

so that the integrand in ( |40| ) is negative over the whole range of integration. 
Because / is continuous we conclude that there exists at least one solution 
y = x a of the equation f(y) = in the range (37) 
luesti 

f'(y) 



The answer to question (ii) can be found by computing f'(y) , which yields 

W'{(e + y)/(e + l))-W f (y) 



1-3/ 



(44) 
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Because W is a monotonic increasing function, and because of the ordering 

» < < 1. («) 



it follows that f'(y) > for all y in the range (37). Thus the answer to 
question (ii) is also positive. 

In conclusion, an annular sector can always be closed into a tube with 
an accompanying axial stretch if its strain-energy function W(x) is strictly 



convex and its thickness is small enough, as prescribed by (38). 



3.1 Thin- walled tubes 



We have shown in the previous section that if the strain-energy function is 
strictly convex, then a unique residually-stressed configuration exists for all 
such energy functions and for any given opening angle and axial stretch. The 
dimensions of this residually-stressed tube depend on the form of the strain- 
energy function, the initial geometry of the unstressed sector, and the value 
of the axial stretch. For thin tubes, some general statements can be made 
about the form of the residually-stressed configuration, as follows. 



Expanding f(y) in (39) as a Maclaurin series in the thickness parameter 

0, and retaining only the leading- 



e, substituting into the equation f(x a ) 
order term, yields the equation 



W (x, 



(46) 



which, by (27), has the unique solution x a = x m = 1/k. Thus, thin tubes 



have the nice property that the residually-stressed configuration B r corre- 
sponds to the minimum of the strain energy. 



Moving on to thicker tubes, we expand equation (36) to the next order 
in e and this yields 



W'{x a )+ l -e (I -,-„)ir"(.r,,)-Tr'(.r„) 



0. 



(47) 



Then we look for x a in the neighborhood of x m , setting x a = x m [1 + e5 + 0(e 2 )}, 
say. Substitution of this into the equation above yields 



Xm5 + \{ l -\ 



W"(x m ) = 0, 



(48) 
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but W" > 0, so that x m 5 = —(k - l)/(2fc), and hence 

x a = x rn [I - e{k - l)/2 + 0{e 2 )] , (49) 

showing that a thicker initial sector leads to a smaller inner radius of the 
residually-stressed tube (a result independent of the particular form of strain- 
energy function within the considered class). 



3.2 Examples 

We now return to the case of tubes with finite thickness. In the preceding 
analysis, e = (B/A) 2 — 1 was used as the thickness parameter. It gives a 
measure of the thickness of the unloaded, unstressed annular sector Bq. It fol- 
lows easily from simple kinematics that the thickness ratio of the residually- 
stressed tube in configuration B r , i say, is related to e by 



- 1 



> £. 



•L a 



(50) 



To make further progress, the strain-energy function has to be specified. 



First, we consider the Mooney-Rivlin model (33). For this material, it is 



found that the defining equation (36) for the residually-stressed state is in- 



dependent of the material constants C and D and does not depend explicitly 
on the axial stretch X z ; it is given by 



In 



e + x a 



- k 2 In (e + 1) + 



(51) 



or, equivalently, in terms of e, 



In (e + 1) - k 2 In (1 + ex a ) + 



e + 1 



1 a 



0. 



(52) 



The constants in these equations now need to be specified. |Matsumoto and 
Sato (2002) measured the inner and outer radii of five tubular segments of 
bovine thoracic aortas, with mean values of a = 11.48 mm and b = 17.45 
mm, respectively, yielding e = 1.310. On cutting these segments longitu- 
dinally, they found an average opening angle of a = 139°, which yields 



k = 1.629. Substituting these values into (52) and solving the equation 
numerically yields x a = 0.5157. This value and (50)2 then yield e = 0.6759 
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Figure 2: Variation of the non-dimensional radial stress a rr and hoop stress 
oqq through an arterial wall modeled as a Mooney-Rivlin material. The 
tube is in a residually-stressed configuration. The opening angle is taken as 
a = 139°. 



which, incidentally, is indeed smaller than l/(k — 1) = 1.590, in line with the 
requirement (38). Finally, the outer value of x is determined as Xb = 0.7110 
by (36)2, and the radial stress reaches its minimum at x = x m = 0.6139. In 
Figure |2j we plot the variations of the resulting (residual) stresses a rr and 
<7gg, divided by the quantity {C + D\ 2 Z )/ (k\ z ), as a function of x. These 
figures are entirely consistent with those produced by other authors, as in, 



for example, Ogden (2003), Raghavan et al. (2004) and Olsson et al. (2006) 



(1997) 



For our second example, we use data from an article by Delfino et al. 
who consider measurements made on human carotid arteries. In 



their experimental, theoretical, and numerical investigation of the mechan- 
ical behavior of arteries, Delfino et al. (1997) model the arterial wall as 



a homogeneous, hyperelastic, isotropic solid, exhibiting residual stress and 
strain-stiffening effects. They use the opening angle method to account for 
the former and the Fung model for the latter. They report the following 
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4.0 mm, and a = 100°. These give here: 



typical values: a = 3.1 mm, b 
e = 0.6649 and k = 1.385. 
surements to the predictions of the Fung strain energy density (35) and find 
the following material parameters: fi = 44.2 kPa, c = 8.35. It follows that 



Delfino et al. (1997) go on to fit pressure mea- 



x a can be determined by solving (36), which reads here as 



Xb W' F {x), 

dx, 

1 — x 



x b 



(e + l)x a 
ex a + l 



(53) 



Numerically, we find that x a = 0.6701. We may then deduce the expected 
value for the radius of the opened sector to be A = a^k\ z /x a = 4.673 mm. 
This comes within less than 5% of the experimentally measured radius, which 



is 4.41 mm (Delfino et al., 1997). 



4 Stability of the residually-stressed configu- 
ration 

In the absence of body forces the equilibrium equation can be written quite 
generally as 

DivS = 0, (54) 

where S is the nominal stress tensor and Div is the divergence operator with 
respect to B$. For an elastic material with strain-energy function W(F) and 
subject to the incompressibility constraint det F = 1, the nominal stress is 
given by 

s-W-^ 1 . («) 

where q is again the Lagrange multiplier associated with the constraint. In 
terms of the Cauchy stress tensor cr the equilibrium equation has the equiv- 
alent form diver = and we note the connection cr = FS. It was the 
specialization of the latter form of the equilibrium equation to the cylindri- 
cal geometry that was used in Section 2.2. In considering the stability of the 
configuration B r it is convenient to work in terms of the incremental form of 



the equilibrium equation (54) and the incremental form of the constitutive 



equation (55), which we now summarize. 
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4.1 Incremental equations 



Let u be a small displacement from the configuration B r and let S be the 
associated increment in S. The (linearized) incremental form of the consti- 



tutive equation (55) is written as 

S = AF + qF -qF-\ 



(56) 



where F and q are the increments in F and q, respectively, and A is the 
fourth-order tensor of fixed-reference elastic moduli, defined by 



d 2 W 



A, 



d 2 W 



" dFdF 1 ~~ am dF ia dF jP 
The incremental form of the equilibrium equation is then 

DivS = 0. 



(57) 



(58) 



Let x denote the position vector of a point in B r . Then, by updating 
the reference configuration from B to B r , we may write the incremental 
equilibrium equation in Eulerian form as 



divSV 



0. 



where 



s = Ar + gr - qi, 



(59) 



(60) 



is the push forward FS of S, F is the displacement gradient gradw, J is the 
identity tensor, and Aq, the push-forward of A, is the tensor of instantaneous 



moduli, whose components are given by A 



Opiqj 



FpaFqflAaiflj . 



For full details 



of these connections, see (Ogden, 1984). The incremental incompressibility 
condition is then 

trr = div-u = 0. (61) 

In the configuration B r it is appropriate to work in terms of cylindrical 
polar coordinates r,8,z and their associated basis vectors e r ,eg,e z . Let 
(u, v, w) be the components of u. Then the matrix of components of T with 
respect to this basis is 

u ,r (w,6» — v)/r U r 

v, r (u + v t e)/r v tZ (62) 
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and the incremental incompressibility condition becomes 



U r H + W 



(63) 



Referred to principal axes of B the non-zero components of «4o are (Og- 



den, 1984): 



"Aonjj X j Xj , 

Aon*, = {KWi - XjWj^/iX* - X 2 j), i ± j, 

*A-0ijji •A-Ojiij -A-Oijij AjWj; i 7^ J: 



(64) 



dW/dXj and W t 



VI 



(no sums on repeated indexes here), where Wj 
d 2 W/dXidXj. 

In the present situation the cylindrical axes are principal axes of B and 
the above apply with % and j running over values r, 9, and z. The compo- 
nents Soij are then easily obtained from (60) and (62). Explicit expressions 
for the incremental equations of equilibrium for a thick-walled tube in cylin- 
drical polar coordinates can be found in Haughton and Ogden (1979). Here, 
however, for purposes of illustration and for simplicity we now specialize the 
incremental equations so that w = and u and v are independent of z, so that 
we consider prismatic incremental deformations. The equilibrium equation 
(59) reduces to two component equations, which can be written as 

(rSorr),r + SW,e — Soeg = 0, [rSo r e),r + Soeefi + SW = 0, (65) 



and from (60) we obtain 



So rr 

Sore 

So6r 

Soee 



I A \ A U + Vfi. 
(Aorrrr + Q) u ,r + Ao rr 60 Q, 



•Aor8r8V,r + \AQ r 6r6 ~ ~rr 
Ua - V 



u , 



A 



Oerdr- 



I A \ U + V fi A 

{Aoeoee + q) h Aon-eeu^ — q, 



(66) 



and the incremental incompressibility equation ( 63 ) reduces to 

u + Va 



0. 



(67) 
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Note that use has been made of the connections 

•Aor06r + q — Ao r g r g — <7 rr = AqOtOt ~ ^ee, (68) 



for which equations (10) and (64) have been utilized. 

Note that the components of incremental traction on a surface r = 
constant are So rr and So r e. 

4.2 Stroh formulation 

We now consider solutions of the form 

{u,v,q} = {U{r),V(r),Q(r)}j n6 , (69) 

where U, V, and Q are functions of r only, so that the solutions are single 
valued, and n = 0, 1, 2, 3, . . . is the circumferential number. It follows that 
the components Sq^j have a similar structure and we therefore write 

£ Qii = Ey(r)e infl , (70) 

say, where Ey, i,j G {r, 8} are functions of r only. Then we obtain 

E rr = (A 0rrrr + p)U' + -A 0rr9e (U + mV) - Q, 

r 

E re = ArOreV + -{Arere - a rr )(inU - V), 

r 

Ti0 r = -Aoer8r(inU — V) + (A()8r6r — ^ee)V , 

r 

E ee = - {Aoeeee + p){U + inV) + A 0rrde U' - Q, (71) 



and the incremental incompressibility condition (67) yields 



U' = --U--V. (72) 



From (71) we obtain 



V = --(1 -a)U + -(l- a)V + (73) 
r r a 

where we have introduced the notations 

a = A 0r 8r8, a = a Tr /a. (74) 
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By using the above and then eliminating T,g r , X^, and Q from the equilibrium 
equations in favor of U, V, £ rr , and H r g, the two equilibrium equations yield 
expressions for (rE rr .)' and (rS^)' in terms of U, V, £ rr , and Y> r g. Then, by 
introducing the four-component displacement-traction vector 77 defined by 



77 = [U, V, irY, rr , ir£ r0 ] 



(75) 



we may cast the governing equations in the Stroh form (Shuvalov, 2003) 

d 



-rj(r) = -G(r)r)(r), 
ar r 

where the matrix G has the Stroh structure, 



(76) 



G 



1 

-n(l - a) 
-i«i2 



-n 
-i(l-a) 

K22 





-1/a 
— i — n(l — cr) 
— n ifl — a) 



(77) 



where 



Kll = 2/3 + 2a(l - cr) + n 2 [7 - a(l - cr) 2 ] , 

K12 = n [2/3 + 7 + a(l - cr 2 )] , 

/? 22 = 7 - a(l - a) 2 + 2n 2 [/3 + a(l - a)], 

and the notations 7 and /3 are defined by 

7 = "Aog r o r , 2/3 = Ao rrrr + Aoggge — 2Ao rr gg — 2Ao r 66r- 



(78) 



(79) 

Note that K n + k 22 = (n 2 + 1)k 12 /?t,. Note also that this form of the Stroh 
formulation is entirely consistent with previous forms obtained in the bending 



of blocks (Destrade et al. , 2009, 2010), or the compression of tubes (Goriely et 



al. , 2008), the only differences being in the actual expressions of the moduli. 



The above coefficients may also be conveniently expressed in terms of 
W(x) through the connections 



a 



2xW'{x) 
k 2 x 2 — 1 



7 = k x a 



(3 = 2x 2 W"(x) + xW'{x) 



a. 



(80) 



and a is given by ( 19 )i and (74). 
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We may now devise a numerical strategy to find out whether a given open 
sector remains stable once closed into a tube, as follows. The prescribed 
geometrical quantities are the measure of the initial sector thickness e, the 
measure of the opening angle k, and the axial stretch X z . The prescribed 
physical quantity is the strain energy density W(x). These quantities are 
required to determine x a and Xb from (36). Then we integrate numerically 



the incremental equations (76), written as 



dx 



77(£) 



2x1 



x) 



-G(x)r}(x), 



(81) 



(using (fl7|)), and we check whether (stability) or not (instability) the inner 
and outer faces of the tube are left free of incremental traction, i.e. whether 
or not 



r){x a ) 



[U(x a ),V(x a ),0,0f 



rt(xb) = [U(xb),V(xb),0,0f 



(82) 



This two-point boundary-value problem may run into numerical stiffness 
problems, especially for thick sectors, and we use the compound matrix 



matrix method (Destrade et al. 2009) to take care of these 



method (Haughton and Orr, 1997 Destrade et al. 2009) and the impedance 



4.3 Example: Mooney-Rivlin tubes 



Motivated by simple experimental observations on sectors made of silicone 
rubbers, see Figure [2j we now focus on tubes made of solids with the Mooney- 
Rivlin strain-energy density (33). We find that 



a 



and 



(CA 



-i 



a 



DX, 



x 
2 



1 

kx' 
'1 - k 2 ) In 



7 = (CX~ + DX z )kx, 2/3 = a + 7, (83) 



x 



x b - 1 



-In 



■t-'b 



1 

+ - - 

x 



1 



(84) 



In the latter equation, x b is determined by solving (51) for x a , a calculation 
which is done independently of C, D, and X z . 



It is now clear upon inspection of the incremental equations (81) that 



once we normalize the third and fourth components of 77 with respect to the 
quantity (CX^+DXg), we end up with a completely non-dimensional system 
of differential equations which is independent of the material parameters C 
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Figure 3: Annular sectors made with rubber silicone; Height: 6 cm; Inner 
radius: A = 3.5 cm; Outer radius: B = 7.2 cm. Two rectangles of acrylic 
glass are glued to the end faces of the sectors in order to bend them. On the 
left: an unstressed open sector with opening angle a = 7r/3 is closed into a 
residually stressed and unstretched tube, and its inner curved face remains 
free of wrinkles. On the right: an open sector with opening angle a = 27r/3 
is bent but cannot be closed without the appearance of axial wrinkles on its 
inner face, one of which is marked by an arrow. 

and D and of the axial stretch X z . Only three parameters are then required 
to determine the stability of closed-up Mooney-Rivlin tubes, namely the 
following geometrical quantities: B/A, the ratio of the initial outer to inner 
radii; a, the opening angle; and n, the mode number. 

Figure [1] shows the variations of the critical opening angle a cr (in degrees) 
with B/A, for different values of n. The lowest curve represents the bifurca- 
tion curve, below which all tubes are stable. In our numerical investigation, 
we found that this corresponds to the n = 9 curve, although the n = 8 and 
n — 10 curves are almost indistinguishable from it. It follows that when 
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a closed-up Mooney-Rivlin tube buckles, it presents at least 8 wrinkles on 
its inner curved face. When the initial sector is thin, B/A — > 0, it can be 
bent into a tube without buckling, so that a CT — > 360° in that limit, for all 
mode numbers. This is consistent with the observation that thin rectangular 
plates with a neo-Hookean strain energy can be bent into full circles without 
encountering bifurcation (Haughton, 1999). 

The radii ratio of the silicone sectors in Figure [3] is B/A = 2.06. We see 
that an opening angle of a = 120° is below the bifurcation curve, so that 
the corresponding sector can be closed into a tube in a stable manner; see 
Point (a) in Figure |4j For an opening angle of a = 240° , we are above the 
bifurcation curve (see Point (b) in Figure [ij), an d wrinkles are thus predicted. 
Finally, the data of Matsumoto and Sato (2002), fitted to the M ooney-Rivlin 

VI + x a i 



strain energy, give B/A 
their measured opening angle of a 



3.2 



In that case, 



= 1.29; see Section 
139° is clearly below the bifurcation 
curve, see Point (c) in Figure [ij and the artery can be "closed" in a stable 
manner. 



5 Conclusion 

We have proved that the opening angle method gives a unique solution to 
the equations of equilibrium provided the strain energy density of the solid 
is a convex function of the geometrical quantity x = 2n\ z ( y r 2 / R 2 )/(2n — a), 
where X z is the axial pre-stretch, a is the opening angle, and r, R are the 
initial and current radial coordinates, respectively. The convexity is indeed 
met by many standard strain energy densities including the Mooney-Rivlin, 
Fung, and Gent models. We have also provided an analysis of incremental 
prismatic deformations superimposed on the finite deformation in order to 
assess the stability of the intact tube formed by the closing of the sector. 

The results are valid for homogeneous isotropic solids and are thus highly 
relevant to the mechanical modeling of the closing of annular sectors into 
tubes made of rubber-like materials. This correlation was confirmed by our 
toy experiments with sectors made of soft silicone. With respect to biome- 
chanics, where the opening angle method is a tool widely used for the esti- 
mation of residual stresses in arteries, we note that there are many articles 
in the literature where arterial walls are modeled as being homogeneous and 



isotropic, as, for example, in Halperin et al. (1991) or Delfino et al. (1997). 



However, it is now well established and accepted that arteries are inhomo- 
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Figure 4: Instability of a closed tube with Mooney-Rivlin strain energy: 
critical opening angle a CT versus initial radii ratio B/A, for different mode 
numbers n. The lowest curve -the bifurcation curve- corresponds to modes 
n = 8, 9, 10 (almost indistinguishable one from another). The curves of 
all other modes are situated above this (n = 2, 3 are shown). The annular 



silicone sector on the left in Figure [3] and the opened artery of Section 3.2 
correspond to stable geometries (Points (a) and (c)). The annular silicone 
sector on the right in Figure [3] gives a point above the bifurcation curve, 
Point (b). 

geneous and anisotropic solids, due, in particular, to a layered structure and 
the presence of collagen fiber bundles. The impact of such inhomogeneity 
and anisotropy on our results remains to be assessed. 
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